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SlilMMARY 


I 

Acoustic  Doppler  Current  Profilers  are  presently  employed  by  the  U  S.  Army  in  a  number  of 
projects.  However,  as  with  most  new  techniques  there  have  been  several  problems  both  with  the 

j 

collection  and  the  subsequent  analysis  of  data.  The  main  objective  of  this  work  is  to  improve  the 
accuracy  of  Acoustic  Doppler  Current  Profilers  by  examining  the  existing  survey  techniques  and 
by  developing  computer  software  to  process  survey  data.  Two  FORTRAN  programs  have  been 
develope\d  during  the  contract  period.  One  program  is  designed  to  combine  data  from  an  external 
differential  global  positioning  system  with  data  from  an  ADCP  to  give  a  more  reliable  discharge 
measurement,  and  the  second  program  is  designed  to  in-fill  gaps  in  records  due  to  missing  ADCP 
^  velociry  data.  In  order  to  test  and  develop  the  software,  data  was  collected  on  the  Columbia 

'  River,  Washington.  This  document  details  the  resulting  software  designed  and  briefly  reports  on 

the  fieldwork  undertaken. 
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I.  CONVERSION  FACTORS  FOR  MEASUREMENT  UNITS 


Because  of  the  parties  involved  in  this  research  project,  namely  The  University  of  Nottingham, 
UK,  and  the  United  States  Government,  a  table  of  conversion  factors  for  Imperial  and  SI  units 
is  given  below: 


FROM 

TO 

MULTIPLY  BY 

Inches  (in) 

Millimetres  (mm) 

25.4 

Feet  (ft) 

Metres  (m) 

0,3048 

Square  Feet  (ft^) 

Square  Metres  (m^) 

0.0929 

Cubic  Feet  (ft^) 

Cubic  Metres  (nT) 

0.02832 

Yards  fyd) 

Metres  (ni) 

0.9144 

Miles  (Mi) 

Kilometres  (km) 

1 .609 

Millimetres  (mm) 

Inches  (in) 

0.0394 

Metres  (m) 

Feet  (ft) 

3.281 

Square  Metres  (m^I 

Square  Feet  (ff) 

10.7639 

Cubic  Metres  (m^) 

Cubic  beet  (fT) 

35.3147 

Metres  (m) 

Yards  (yd) 

1.0936 

Kilometres  (km) 

Miles  (Mi) 

0.621 
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1.  INTRODUCTION 


Acoustic  Doppler  Current  Profilers  (ADCPs)  are  "state-of-the-art"  instruments  that 
measure  water  velocities  acoustically,  in  particular  the  Doppler  effect.  I'he  Doppler  effect  is  the 
change  in  observed  sound  pitch  that  results  from  relative  motion  between  a  source  and  observer. 
The  major  advantages  of  ADCPs  over  more  traditional  methods  of  river  survey  are  their  speed 
of  use,  quality  of  data,  resolution  of  measurement  and  their  ability  to  be  used  onboard  moving 
vessels  to  measure  river  discharges  (Didden,  1 989) 

However,  due  to  their  nature  of  operation  ADCPs  produce  large  quantities  of  data  that 
requires  extensive  post-survey  processing  in  order  to  be  converted  into  a  usable  format. 
Additionally  when  being  used  onboard  a  moving  vessel  and  under  certain  conditions  (i  e.  high 
flows)  their  ability  to  "bottom-track"  breaks  down  (Miinchow  et  al.,  1995).  The  ADCP  records 
the  velocity  of  the  channel  bed  as  well  as  the  water  velocity  and  makes  a  correction  for  vessel 
motion  during  measurement.  In  order  to  overcome  this  "bottom-track"  problem  new  survey 
techniques  and  post-survey  data  processing  software  have  been  designed.  The  "bottom-track" 
facility  is  replaced  by  an  external  Differential  Global  Positioning  System  (DGPS)  to  monitor  and 
record  vessel  motion. 

During  this  study  test  data  sets  were  collected  using  an  ADCP  in  conjunction  with  a 
DGPS  to  measure  water  velocities  and  river  discharge  on  the  Columbia  River,  upstream  and 
downstream  of  the  Bonneville  Dam,  Washington.  GPSQ  a  DOS  based  FORTRAN  computer 
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recorded  from  the  DGPS. 

In  addition  to  the  "bottom-track"  problem,  ADCP  data  can  be  prone  to  missing  ensembles 
in  a  given  cross-section,  f'hese  missing  data  give  rise  to  accumulated  errors  in  discharge  and 
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positioning  data  when  the  ADCP  is  being  used  in  a  standard  mode  (i.e.  without  an  external 
DGPS),  During  this  study  a  second  DOS  based  FORTRAN  program  (ADCPFILL)  was  developed 
to  fill-in  missing  data. 


2.  PROJECT  CHRONOGRAPHY 


2.1  ITINERARY 
1®^  September  1995: 

Mr  Richaidson  and  Prof  Thome  ai rived  in  Vicksburg,  Mississippi.  After  initial  meetings 
Mr  Richardson  began  work  at  Waterways  and  met  with  Mike  Trawle,  Tirn  Fagerburg  and 
Thad  Piatt  to  outline  possible  areas  of  concern  with  present  ADCP  data  collection 
techniques  and  analysis. 

2"*  September  to  30‘*'  Septemlier  1995; 

Bottom-track  enors  were  highlighted  as  a  major  problem  in  the  ADCP  river  gauging 
technique,  and  subsequently  Mr  Richardson  began  developing  a  method  for  calculating 
discharge  using  a  Differential  Global  Positioning  System  to  replace  bottom-track.  He  also 
started  designing  general  code  for  GPSQ  software. 

1“  October  to  14“‘  October  1995: 

Field  work  in  Washington,  on  Columbia  River.  Mr  Richardson  and  Thad  Pratt  collected 
sample  data  sets  necessary  for  software  development  and  refined  survey  methods  for 
collecting  and  logging  DGPS  data. 

16**'  October  to  I"*  November  1995; 

Mr  Richardson  and  Thad  Pratt  returned  to  Vicksburg  to  refine  and  complete  software 
development  for  GPSQ  Problems  of  missing  ADCJP  data  in  all  ADCP  applications  became 
apparent  -  ADCPFILL  was  designed  to  overcome  this 
2'"*  November  1995: 

Mr  Richardson  returned  to  U.K. 
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2.2  DATA  (OM.l  ( HON 

111  order  to  ovcicoiiie  the  inlicrcnt  notlonv-track  problem  Aithin  the  AIX'P  setup,  an 
alternative  to  bottotii-lracking  must  be  utilised  to  monitor  boat  motion  An  ideal  alternative  is  a 
Differential  (ilobal  Positioninti  System  (DGPS).  A  DGPS  aehieves  reasonable  accuracy  in 
absolute  position  and  it  is  also  capable  of  greater  precision  when  monitoring  relative  positions  or 
velocity  (PAP24',  l‘)'-)1) 

A  trial  data  set  was  collected  on  the  Colutnbia  River,  Washington  upstream  and 
downstream  of  the  Bonneville  Bam  Data  was  collected  using  an  RDl  Instmments  ADCP  with 
an  acoustic  frequency  of  600  kBz  in  conjunction  with  a  Trimble  Navigations  DGPS,  A  test  data 
set  was  needed  for  a  number  of  reasons.  First,  it  was  necessary  in  order  to  develop  the  software 
required  for  the  subsequent  corrections  applied  to  the  raw  ADCP  data.  Second,  the  technique  of 
using  a  DGPS  in  conjunction  with  an  ADCP  for  gauging  discharge  required  verification,  funally, 
a  survey  method  for  setting  up  and  logging  data  from  the  DGPS  was  required. 


Plate  1 .  Test  Site  on  Columbia  River,  Washington 
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Figure  1  Location  Map  of  Test  Data  Site,  Bonneville  Dam,  Columbia  River,  Washington. 

2.3  RESULTS 

Figure  1,2  shows  results  obtained  at  one  cross-section  dur.ng  the  test  survey  downstream 
of  Bonneville  Dam.  Figure  L2(a)  shows  results  from  the  ADCP  used  in  bottom-track  mode,  with 
a  discharge  calculated  of  130143  cfs.  Figure  1 .2(b)  shows  results  obtained  using  the  ADCP  in 
conjunction  with  the  DGPS,  with  a  discharge  calculated  using  'GPSQ'  of  191750  cfs.  It  is  obvious 
from  the  plots  of  measured  velocities  in  figure  1 .2  that  the  velocity  field  shown  is  significantly 
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different  dependent  on  the  velocitv  reference  used  Also,  the  discharge  calculated  using  the  IDOPS 
as  a  velocity  reference  is  47%  greater  than  that  when  using  bottom-track.  The  'bottom-track 
declination'  calculated  by  'GPSQ'  at  the  time  of  data  processing  is  21  degrees  (the  bottom-track 
declination  is  the  difference  between  the  survey  line  direction  calculated  using  bottom  track  data 
and  the  survey  line  direction  calculated  using  the  DGPS  data). 

The  discharge  results  may  be  explained  in  part  by  the  theory  that  bottom-track  errors  lead 
to  an  under-estimation  of  the  water  velocities  and  therefore  discharge,  through  an  under¬ 
estimation  of  vessel  velocity.  The  bottom-track  declination  of  21  degrees  may  be  partly  due  to 
very  localised  magnetic  affects  on  board  the  survey  vessel  interfering  with  the  interna!  ADCP 
compass.  This  affects  the  measurement  of  velocities  in  two  ways.  First,  the  bottom  track 
measurement  of  the  vessel's  speed  and  direction  is  affected.  Second,  the  measurement  of  water 
velocities  made  relative  to  the  ADCP  is  affected.  When  these  two  measurements  are  combined 
to  calculate  an  absolute  water  velocity,  the  affect  o^^"  local  magnetic  declination  is  compounded. 
The  DGPS  is  not  susceptible  to  the  affects  of  local  magnetic  declination  and  can  therefore 
eliminate  the  affect  from  the  measurement  of  vessel  velocity. 


distance  along  transect  (ft) 

discharge:  130143  cfs 


0  100  200  300  400  500  600  700  800  900 


distance  along  transect  (ft) 
discharge:  191750  (cfs) 

Figure  1.2  Measured  velocities  (ft/s);  (a)  'Bottom-track'  referenced  and  (b)  GPS  referenced 
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3.  GFSQ 


3.1  INTRODUCTION 

When  ADCPs  are  used  on  board  moving  vessels  for  river  discharge  measurements,  they 
are  generally  configured  to  record  the  velocity  of  the  bed  (bottom  track)  and  the  velocity  of  the 
water  relative  to  the  ADCP.  The  ADCP  then  combines  the  two  velocities  to  give  an  absolute 
water  velocity.  Under  certain  conditions  ADCP  bottom  track  data  can  become  unreliable  and  thus 
introduce  errors  into  the  measurement  of  water  velocities  and  discharge.  ADCP  bottom  track  is 
limited  by  the  convex  arrangement  of  the  ADCP  transducer  heads,  making  it  impossible  for  the 
ADCP  to  measure  the  entir  e  water  column.  When  the  ADCP  is  measuring  what  it  assumes  to  be 
the  bottom,  it  is  in  fact  measuring  a  moving  surface  of  wat  /sediment  just  above  the  bed.  Under 
normal  conditions  the  echo  fiom  the  hard  bottom  surface  is  enough  to  overcome  this  problem, 
but  under  higher  flow  conditions,  when  sediment  is  being  entrained  into  the  water  and  the  bed  is 
mobile,  errors  are  introduced.  Generally,  bottom  track  errors  of  this  nature  will  lead  to  an  under¬ 
estimation  of  the  discharge  due  to  the  ADCP  using  a  velocity  reference  that  is  moving  in  the 
downstream  direction 

'GPSQ'  is  a  DOS  based  program  written  in  FORTRAN  whicli,  together  with  an  amended 
survey  technique,  will  re-calculate  ADCP  measured  velocities  and  discharge  It  does  this  by 
utilising  data  from  an  external  differential  global  positioning  system  (DGl^S),  to  replace  bottom 
track  and  measure  boat  motion.  GPSQ  uses  ADCP  velocity  data  that  is  taken  from  the  RDI- 
TRANSECT  program,  in  the  usual  ASCII  processed  data  format,  with  the  velocity  reference  set 
to  "none".  IXfPS  files  are  then  combined  with  the  ADCP  ASCII  files  to  correct  for  boat  motion 
and  produce  correct  velocities  and  discharge. 
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3.2  METHOD 


GPSQ  works  in  three  modules.  The  program  modules  are: 

( 1 )  DGPS  data  processing  module 

(2)  Velocity  data  correcting  mouule 

(.3)  Discharge  calculating  module 

3,2.1  DGPS  Data  Processing 

The  tirst  module  of  GPSQ  reads  the  GPS  navigation  tile  and  calculates  the  boat  motion 
and  position  for  each  ADCP  ensemble.  The  module  will  read  GPS  data  in  the  following  format; 


I'.NSKMUU':  I0.rs  46K524ii 

10:01:56  102465X.6I  7IXK68.14 
10:01 :57  1624656.07  71X866.96 
10:01:58  1624654.36  718X66.70 
10:01:59  1624652.71  71X867.11 

liNSlCMlll.K  10:»9  It' riMl' 4685605 
1 0:02:00  1 62465 1 .46  7 1 8868. 1 4 
10:02:01  1624650  64  718869.64 
1(1:02:02  1624649.27  71X871.14 
10  02:03  1624648.40  718873.53 
10:02  04  1624648.40  71X873.53 

l'.NSi:.MIII.l.  1040  I’CriMi: 46X6034 
10:02:05  1624(.47.33  71X876.20 
10:02:06  1624646.06  718878.92 
10:02.07  16'i4(,44.58  71X881.78 
10:02:08  1624642.61  71X885.42 


The  ensemble  header  line  gives  an  ensemble  number  followed  by  an  ensemble  start  time 
(PCTIME)  in  hundredths  of  a  second,  The  GPS  data  lines  give  a  position  time 
(hours:minutes:secoEids)  followed  by  an  Easting  and  Northing  The  GPS  data  processing  module 
reads  this  file  into  an  array  of  ensemble  number  with  PCTIME  from  the  end  of  the  ensemble  to 
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the  start  of  the  transect,  and  an  array  of  GPS  East  and  North  with  time  from  the  start  of  the  first 
ensemble.  The  two  arrays  are  then  inalched  in  terms  of  time  to  give  each  ensemble  number  a 
position  at  the  end  of  the  ensemble  and  an  average  boat  velocity  duiing  the  ensemble  (a  velocity 
is  calculated  for  each  GPS  record  and  then  averaged  over  the  ensemble  period). 

The  next  segment  of  this  module  then  caJculates  a  survey  line  definition  based  on  the  GPS 
data.  The  start  and  end  points  of  the  GPS  record  are  used  as  the  start  and  end  points  for  the 
survey  line.  A  definition  of  the  sutvey  line  is  necessary  for  the  discharge  calculation  method 
described  below.  A  distance  along  the  survey  line  (course  made  good;  CMG)  and  a  width  is 
calculated  for  each  ensemble,  together  with  a  corrected  position.  GPS  positions  are  then 
translated  to  survey  line  positions  using  the  following  equations: 


E 

= 

(m  n c  -  c.m)  /  (m'^  ^  1)  . 

. (Eqn3,l) 

N 

m.  E  t  -  c  . 

. (Eqn3,2) 

where: 

E 

Corrected  liasting 

N 

Corrected  Northing 

c 

(rPS  Easting 

n 

(IPS  Northing 

m 

Survey  line  gradient  {dist.  North  /  dist.  East) 

c 

Survey  line  intercept 

The  output  of  the  (iPS  processing  module  is  in  the  format  shown  below.  Velocities  and  course 
made  good  data  are  given  in  the  selected  data  output  units,  US  customary  or  SI.  The  GPS 
processed  data  file  is  a  temporaiy  file  which  is  removed  at  the  end  of  the  program  execution. 
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I'NS 

CNKr 

I'ASTINd 

NORTIIINCJ 

!■:  VI- 1. 

N 

WIDrfl 

I0.1S 

2.96 

1624658.00 

71K86K.90 

.S.Kl 

-5  4.1 

1.44 

10.^9 

18.76 

1 624656.00 

718871.20 

2.5. 9.1 

-.1.69 

5  67 

1040 

52.84 

1624651.00 

71X882.40 

42.18 

-5.20 

9.05 

1041 

9S.7:i 

162464.S.0O 

718894.80 

49.46 

-6.25 

9.54 

1042 

155.82 

1 6246.16  00 

7189I0.U/ 

64.82 

-745 

9.78 

104.1 

206..U 

1 624629.00 

718924.70 

46.06 

-7.48 

9,68 

1044 

252.05 

162462.4.00 

7 1X9.16.00 

55.45 

-5.97 

10.25 

1045 

.in. 51 

1624614.00 

718952.60 

67.55 

-5,54 

11.29 

3.2.2  Velocity  Correction 

The  next  module  in  (iPSQ  takes  the  processed  GPS  data  and  combines  it  with  the 
processed  ADCP  ASCII  data,  to  give  corrected  water  velocities.  The  ADCP  processed  ASCII 
tiles  are  in  the  format  shown  in  Appendix  A,  ^ 

For  each  measurement  bin  the  measured  velocity  direction  is  corrected  for  the  local 
magnetic  declination  (P)  as  input  by  the  user. 


P  -  P-'l'  . 

v/here: 

P  Local  magnetic  declination 

p  fiuc  Al)(d’  direction 

i|i  Measuicd  AnC.’P  direction 


(riqn3,3) 


The  non-rd'crenced  East  and  North  ADC!P  velocities  and  then  corrected  for  boat  motion  to  give 
an  absolute  watc!  velocity. 


V,.  '  V^„,, 


(Eqn  3,4) 


where; 


"  Absolute  water  velocity 
V„  Boat  velocity  given  by  DGPS 

'^Aix  p  "  ADCP  non-referenceH  measured  water  velocity 

A  depth  is  calculated  for  each  ensemble  using  either  ADCP  beam  depth  data  given  in  the 
ADCP  ASCII  file,  or  external  fathometer  depth  data  in  a  format  as  yet  untested  in  the  field.  The 
ensemble  depth  is  used  to  calculate  an  ADCP  range.  Velocity  data  outside  of  this  range  is 
discai'ded.  The  range  of  acceptable  ADCP  data  is  given  by  the  ollowing  equation: 


R 

where: 

R 

0 


(O.,  -  I^Aix  p)  .  cose  . (EqnlTS) 

ADCP  range 
Ensemble  depth 

ADCP  transducer  depth  (given  in  ADCP  ASCII  file) 

ADCP  transducer  beam  angle  (INPUT  by  user) 


fhe  output  from  the  velocity  data  conection  module  is  given  in  the  file  VHl.OCn’Y.DAT 


which  IS  created  in  the  working  directoi^.  Vl'.l  XX’fl'Y.DAT  is  in  the  following  format:  • 


id,  «)K, 

disl. 

cast, 

iiiRtli, 

vvltWi, 

UlvUi, 

liin, 

IIIBfi, 

dim. 

cart, 

nt»rth, 

vert. 

errur 

,  (i!)  avji  ndq)  raoj'c 

1  am 

2.97 

1024038. OU 

718868.90 

.5.81 

24.92 

5.00 

5.61 

297.30 

-4.99 

2.57 

.00 

1.70 

247.25 

'*3. 48 

1  nm 

2.97 

1024^>5«  00 

718868.90 

5.81 

24.92 

6.70 

5.49 

288.80 

-5.20 

1.77 

-.10 

-.60 

2.35,50 

23.48 

I  I03X 

2.97 

ir»2465«00 

718868.90 

.5.81 

2492 

K..3(} 

5.19 

282.71 

-5.06 

1.14 

.00 

-.40 

206.25 

23.48 

1  I03K 

2.97 

17324658.00 

718868,90 

5  81 

24  92 

9.90 

5.37 

277.X7 

-5.32 

.74 

20 

-.60 

180,25 

2.1.48 

t  in3K 

2.97 

1  (>24658  0(1 

718868.90 

'81 

24.92 

1 1 .60 

4.52 

282.69 

-4.41 

.99 

.10 

-.10 

172.50 

23.48 

1  I03X 

2.97 

I(>24658.0() 

718868,90 

5.81 

24.92 

13.20 

5.25 

291.11 

-4.90 

I.H9 

.20 

-.20 

16,3.00 

23.48 

1  103K 

2.97 

1624658.(10 

718868  90 

5.81 

24.92 

14.90 

5  46 

290.98 

-5.10 

1.95 

.20 

-.70 

159.25 

23.48 

]  103X 

2.97 

1 624658. (X) 

718X68.90 

3.81 

24.92 

16.50 

5.65 

297.86 

-4.99 

2.64 

.10 

-.90 

153.51) 

23.48 

I  HBK 

2  97 

1 624658. (K1 

718868.90 

5.81 

24.92 

IK.IO 

5.86 

303.79 

-4.87 

.3.26 

-.20 

.00 

150.50 

23,48 

I  1038 

2.97 

162465H.(M) 

/18868.90 

5.81 

24.92 

19.80 

5.14 

289.47 

-4.85 

1.71 

-  10 

-  10 

165.50 

2.3.48 

2  1039 

18.77 

16247»56.00 

71X87.3,20 

25.94 

25  50 

5.00 

5.25 

.307.51 

-4.17 

3.20 

,20 

-.10 

21 1. 00 

24.02 

2  1039 

IH.77 

1624656.00 

71887.3,20 

25.94 

25.50 

6.70 

5.42 

304.90 

-4.45 

.3.10 

.10 

,10 

204.75 

24,02 
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Any  bad  bins  of  measurement  within  the  ADCP  range  arc  ignored  and  are  not  recorded 
in  the  corrected  velocity  file. 


3,2.3  Discharge  Calculation 

A  discharge  is  calculated  based  on  the  survey  line  defined  by  the  GPS  data.  Each  ensemble 
is  assigned  a  width  and  position  along  the  survey  line  and  a  discharge  is  calculated  through  the 
ensemble.  For  each  ensemble  a  top  layei  discharge  (estimated),  middle  layer  discharge  (measured) 
and  a  bottom  layer  discharge  (estimated)  are  calculated.  Two  options  arc  available  for  the  top 
layer  discharge  estimation,  the  power  or  con.stant  method,  fhe  missing  bottom  layer  discharge 
is  calculated  using  the  power  method 


TOP  LAYER  DISCHARGE  ESTIMAI  ION  USING  THE  POWER  METHOD 

The  power  law  equation  for  describing  velocity  profiles  in  open  channels  (Chen,  1991 : 
Simpson  «&  Oltmann,  1990)  is  given  as: 


where: 


9.5.(Z/Z„y’  . (Eqn3.6) 


Velocity  at  height  Z  alxwc  fhe  bed 


Shear  velocity 

Height  above  channel  bed 

Bottom  roughness  height 


Power  Law  exponent 


which  can  be  re-arranged  to  give: 
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(9,5  .  U^/Z,”)  .  Z 
(9.5  .U*/Z„'')  .. 


(Eqn  3.7) 
(Ec|n3.8) 


where  a'  represents  the  non-ADCP  parameters,  then 


(Eqn  3.9) 


For  each  measurement  bin,  GPSQ  solves  equation  3.9  to  give  a  value  of  a'.  Within  each  ensemble 
the  a'  values  are  averaged  (a'  avg),  and  the  average  is  then  used  to  solve  for  the  missing  top  layer 
velocity.  The  top  layer  discharge  is  then  given  by. 


where: 


a’  avg  .  ZV  W  ,  D,, 


Top  layer  discharge 


Corrected  ensemble  width 


Blank  depth 


(Eqn  3.10) 


TOP  LAYER  DISCHARGE  ESTIMA  TION  USING  THE  CONSTAN  T  MI3THOI) 

The  constant  method  for  estimating  the  missing  dischaigc  near  the  surface  is  the  simpler 
of  the  two  techniques.  The  constant  method  docs  not  follow  accepted  hydraulic  principles  of 
velocity  profile  description,  however,  under  certain  circumstances  it  seems  to  give  a  more 
accurate  picture  of  near  surface  properties.  The  niethou  simply  takes  the  velocity  from  the  first 
good  measurement  bin  in  the  profile  (V^,,,)  and  extrapolates  to  the  suiface. 


W  .  . (Eqn  3.11) 
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MIDDLE  LAYER  MEASURED  DISCHARGE 

The  discharge  through  each  measurement  bin  is  given  by  the  following  equation: 


Qbin 

where: 

Qhu, 

V 

a 

P 

H 


V  .  cos(a  -  p;  .  W  .  H . (Eqn  3.12) 

Bin  discharge 

Corrected  velocity  magnitude 
Survey  line  direction 
Corrected  velocity  direction 
Bin  size 


The  bin  size  (H)  is  equal  to  the  distance  between  the  bin  being  considered  and  the  last  good 
measurement  bin,  This  is  usually  equal  to  the  configured  bin  size  but  will  change  if  any  bins  are 
missing  due  to  bad  data.  I’he  middle  layer  discharge  is  then  equal  to  the  sum  of  ail  measurement 
bin  discharges, 


BO  I'J’OM  LAYER  EST1MA3  ED  DISCHARGE 

The  missing  dischai  gc  outside  of  the  ADCP  range  at  the  bottom  of  the  velocity  profile  is 
calculated  using  the  power  law  equation  (Eqn  2.5).  The  bottom  layer  discharge  is  given  by: 


blm 


where: 


Bottom  layer  discharge 


ri^np  -5  n;\ 


O 
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NEAR  SHORE  DISCHARGE  ESTIMATION 


GPSQ  can  be  used  to  estimate  any  discharge  that  is  not  measured  close  to  the  banks  at 
the  beginning  and  end  of  each  transect.  GPSQ  uses  the  same  ratio-interpolation  method  employed 
by  TRANSECT  to  estimate  near  shore  discharge.  The  method  uses  the  first  and  last  known 
velocities  and  assumes  a  tri-angular  area  between  the  bank  and  transect.  The  missing  near  shore 
discharges  are  given  by; 


Q 


Rhcrc 


where: 


Qnluiic 

V... 


L 


0,707  .  V,„  .L  .D„./2 


(Eqn3.14) 


Near  shore  discharge  (start  or  finish) 
Nearest  ensemble  depth  averaged  velocity 
Distance  to  bank 


D„,  =  Nearest  ensemble  depth 


rhe  total  discharge  lor  the  transect  is  then  given  by  the  sum  of  the  top  layer  discharges, 
measured  dischaiges,  bottom  layer  dischaiges  and  the  near  shore  discharges.  The  discharge  output 
file,  whose  name  and  location  is  specified  by  the  user,  is  in  the  following  format: 

lo^i  lay‘Tdisdi.irgccsli!natii)ii  iiKtluxl:  power 
!  piiwtT  oo-ciliaait  1 .6670001 '.-Ol 


ctisatildc,  dist.  easting 

notiliinp. 

widlli, 

(kpth,  liin  tcif)  q,  Umri.  total  q 

uus 

2.97 

Ki2465«.0() 

7imx.9() 

5,KI 

24.92  -3.40.10  -95.00  -1K.51  -150,21  -2.34 

!(U‘; 

IS. 77 

1  f)24656.00 

7IHX7:?.20 

2^.94 

25,50  -9.39.54  .203.17  -52.39  -1255.10  -1 

.4.1 

1()40 

52.S<f 

16240.^1.00 

7188H2  40 

42.  IX 

30.50  -2303. IH  -547.49  -130.89  -3041.50  - 

1.7X 

1041 

9H.74 

1024045.00 

7IK«‘M.K0 

49.46 

37  93  4173.10  088.60  204.95  -5000.70  - 

1.81 

1042 

155.S.1 

10240:10.00 

71X7 1  (7  10 

64. X2 

44.30  -7/50.80  -1012.12  -354.03  -9117.00 

-2.00 

Uilal  diseJiarj^c  -18^21  VM)(dN) 
tcilfll  tTtm-stxlioo-Tl  aioji  n'(07() 


4.  ADCPFSLL 


4.1  INTRODUCTION 

ADCPFILL  is  a  DOS  based  program  which  will  read  processed  ADCP  files  from  the  RDI 
TRANSECT  program  and  write  over  bad  ensemble  data  Bad  ensemble  data  will  occur  when  two 
or  more  ADCP  transducer  beams  do  not  return  good  data  for  any  leason  Bad  ensemble  data 
affects  the  overall  calculation  of  discharge  and  the  ability  of  the  ADCP  to  calculate  the  position 
of  each  ensemble,  in  terms  of  distance  travelled,  displacement  east  and  north  and  straight  line 
distance  to  starting  point  (ami  se  made  good).  Therefore,  after  one  bad  ensemble  the  bottom  track 
and  discharge  data  for  every  subsequent  ensemble  has  a  built-in  error,  which  accumulates  after 
every  bad  data  segment. 

4.2  METHOD 

During  playback  of  raw  ADCP  files,  an  ASCII  file  can  be  produced  (using  the  ALT-A 
command)  containing  position  and  velocity/discharge  data  for  each  ensemble.  This  ASCII  file 
contains  a  header  followed  by  velocity  data  for  each  ensemble  in  the  transect.  The  ensemble 
header  contains  the  time  at  the  start  of  the  ensemble,  bottom  track  data,  depth  data  for  each 
transducer  beam,  discharge  data  and  various  other  survey  parameters  (see  appendix  A).  The 
discharge  data  contained  in  the  ensemble  header  can  be  used  as  a  test  for  bad  segments  of 
ensembles.  Bad  ensemliles  show  no  increment  in  the  measured  discharge,  whereas  good  ensembles 
show  an  increase  in  the  cumulative  discharge  measurement. 

Bin  data  for  the  ensemble  is  then  given  in  the  form  of  depth,  velocity  magnitude,  velocity 
direction,  east  velocity,  north  velocity,  vertical  velocity,  error  velocity,  echo  intensities, 
percentage  good  data  (based  on  the  number  of  beams  used  in  the  solution  of  the  equation  of 
motion  of  the  water)  and  bin  discharge.  For  bad  bins  of  measurement,  data  values  of -32768  are 
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attributed  to  velocities  and  2147483647  are  attributed  to  discharge.  Bad  measurement  bins  are 
the  result  of  either  the  bin  being  outside  the  ADCP  range  or  more  than  one  beam  returning  a  bad 
signal  from  that  particular  range  of  depths. 

Once  bad  ensembles  have  been  identified,  ADCPFILL  uses  a  simple  linear  interpolation 
to  re-calculate  ensemble  bottom  track  data,  beam  depths,  pitch,  roll,  coirected  heading,  ADCP 
temperature,  bin  velocity,  echo  intensities  and  percentage  good  data.  Discharge  data  is  not  filled 
by  interpolation,  but  calculated  separately.  Ensemble  time  and  number  are  not  affected  by  bad 
velocity  measurements  and,  therefore,  do  not  need  recalculatiitg.  ADCPFILL  uses  the  following 
equation  to  interpolate  between  good  ensembles: 


ki 

— 

k,(n2-ni)  +  k2(ni-n,)  /  (Oj-n,) . 

where: 

k 

~ 

Parameter  under  consideration 

n 

-- 

Ensemble  number 

i 

= 

Bad  ensemble 

1 

Last  good  ensemble  before  bad  data  segment 

1 

Next  good  ensemble  after  bad  data  segment 

(Eqn  4.1) 


Once  the  bottom  track  and  velocity  data  for  each  bad  ensemble  has  been  corrected,  the  bin  and 
ensemble  discharges  can  then  be  re-calculated.  ADCPFII  ,L  uses  the  same  moving  boat  method 
of  discharge  calculation  employed  by  TRANSECT  (Simpson  and  Oltmann,  1990;  Gordon,  1989). 
fhe  discharge  for  any  given  measurement  bin  is  calculated  as  follows: 


Qw„  ^  X,  .  hi.  At 


(Eqn  4.2) 
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where 


Q,,„,  '=  Fiin  discharge 

Xj  Vector  cross  product 

h,  Bin  height 

At  Time  elapsed  during  ensemble 

The  vector  cross  product  (XJ  is  given  by  the  following  formula: 


■  ^watcr)  '  (^boil  ■  ^w»lcr)  •  •  •  • 

where; 

Eh,„i  --  East  velocity  of  the  survey  vessel 
Nh,,,,  North  velocity  of  the  survey  vessel 

E^„,<,r  =-  East  velocity  of  the  water 

Nvaiet  ^  North  velocity  of  the  water 


(Eqn  4.3) 


The  accumulated  errors  in  measured  dischaige  and  bottom  track  data  are  calculated  for  each  bad 
segment  of  ADCP  data  and  then  added  to  subsequent  ensembles  to  correct  for  accumulated 
errors. 

The  final  output  of  ADCI’FILE  is  in  the  same  format  as  the  original  ADCP  data  file,  and 
as  such  can  be  used  by  other  software  packages  designed  for  ADCP  ASCII  data.  3Te  numbers 
of  any  bad  ensembles  eon  ected  are  displayed  on  the  scieen  at  tlie  time  of  processing. 
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5.  CONCLUSIONS  &  RECOMMENDATIONS 


At  conception,  this  project  was  designed  to  address  the  general  problem  of  ADCP  data 
management  and  analysis  on  the  Mississippi  River,  in  particular  to  outline  the  basis  for  a  usable 
engineeiing  database  of  ADCP  data.  However,  under  certain  flow  conditions  a  more  iiindamental 
problem  was  identified,  that  of  actually  collecting  good  quality  data  using  the  ADCP.  Using 
experience  gained  working  in  Bangladesh  with  ADCPs  on  the  Brahmaputra  River,  a  new  survey 
and  data  processing  technique  was  designed  and  tested. 

In  the  main  body  of  this  report  two  computer  programs  are  detailed,  GPSQ  and 
ADCPFILL,  which  were  developed  during  the  study  to  overcome  two  of  the  basic  problems  of 
ADCP  data  collection,  bottom-track  errors  and  blank  or  missing  data.  These  programs  were 
validated  and  found  to  work  well  with  the  sample  data  collected  on  the  Columbia  Rivet ,  however 
they  require  further  development  in  order  to  be  more  flexible  for  general  use  and  eventually  to  be 
incorporated  into  a  standard  method  of  river  discharge  measurement. 

The  significant  differences  between  ADCP  velocity  and  discharge  data  collected  using 
bottom  tracking  and  DGPS,  would  seem  to  pose  more  questions  then  they  answer,  Additional 
work  is  needed  analysing  the  data  collected  during  the  test  survey,  in  particular  looking  at  the 
affects  of  local  magnetic  declination  and  the  possible  incorporation  of  a  gyro  compass  into  the 
survey  method. 

Further  improvements  to  the  survey  methods  could  be  achieved  through  the  addition  of 
an  external  fathometer  gauge  to  measure  cross-sectional  bathymetry.  This  would  require  further 
soffwaie  development  and  field  testing.  At  present  depth  data  is  obtained  from  the  ADCP  beams, 
a  fathometer  would  give  a  more  detailed  and  accurate  measurement,  thus  improving  confidence 
in  the  discharge  calculation.  Further  work  is  also  needed  on  procedures  to  ensure  correct 
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matching  of  the  ditferential  global  positioning  system  file  and  the  ADCP  file,  in  order  to  ensure 
that  data  taken  from  the  two  files  are  combined  accurately  in  time. 
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APPENDIX  A 


ADCP  ASCII  Output  Format 


Row 


Fi  eld 


Description 


I 


1 

2 

3 

4 

5 

6 

7 

8 
g 

10 

11 

12 

13 


ENSEMBLE  TIME  -  Year  (at  start  of  ensemble) 

-  Month 

-  Day 

-  Hour 

-  Minute 

-  Second 

-  Hundredths  of  seconds 

ENSEMBLE  NUMBER  (or  SEGMENT  NUMBER  for  processed  or  averaged  raw  data) 

NUMBER  OF  ENSEMBLES  IN  SEGMENT  (if  averaging  ON  or  processing  data) 

PITCH  -  Average  for  this  ensemble  (degrees)  , 

ROLL  -  Average  for  this  ensemble  (degrees) 

CORRECTED  HEADING  -  Average  AOCP  heading  +  heading  offset  -  magnetic  variation 
ADCP  TEMPERATURE  -  Average  for  this  ensemble  (*C) 


2 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 


BOTTOM-TRACK  VELOCITY  -  East(+)/West(-) ;  average  for  this  ensemble  (cm/s  or  ft/s) 

-  North(+)/SQuth(-) 

-  Vertical  (upW/downC-]) 

-  Error 

WATER  LAYER  VELOCITY  -  East(+)/West(-) ;  average  for  this  ensemble  (cm/s  or  ft/s) 

-  NQrth(+)/South(-) 

-Vertical  (up[+]/dQwn[-l ) 

-  Error 

DEPTH  READING  -  Beam  1  average  for  this  ensemble  (m  or  ft.  as  set  by  user) 

-  Beam  2 

-  Beam  3 

-  Beam  4 

-  Beam  S  (if  present) 


3  1  TOTAL  ELAPSED  DISTANCE  -  Through  this  ensemble  (from  bottom-track  data;  in  m  or  ft) 

Z  ELAPSED  TIME  -  Through  this  ensemble  (in  seconds) 

3  DISTANCED  TRAVELLED  NORTH  (m  or  ft.  as  set  by  user) 

4  DISTANCED  TRAVELLED  EAST  (m  or  ft.  as  set  by  user) 

5  COURSE  HADE  GOOD  -  Through  this  ensemble  (freen  bottom-track  data) 


4 


S 


6 


1  NAVIGATION  DATA  -  Latitude;  average  for  this  ensemble  (degrees) 

2  -  Longitude;  average  for  this  ensemble  (degrees) 

3  -  North(+)/Sauth(-}  ship  movement;  average  for  this  ensemble  (cm/s  or  ft/s) 

4  -  East(+)/West(-)  ship  movement;  average  for  this  ensemble  (cm/s  or  ft/s) 

5  -  Total  distance  Iravellnd  through  this  ensemble  (on  or  ft) 

1  DISCHARGE  VALUES  -  Kiddle  part  of  profile  (measured);  m’/s  or  ftVs 

2  -  Top  part  of  profile  (estimated);  mVs  or  ft’/s 

3  -  Bottom  part  of  profile  (estimated);  ra‘/s  or  ftVs 

4  -  Start-shore  discharge  estimate;  mVs  or  ftVs 

5  -  Starting  distance  (boat  to  shore);  in  or  ft 

6  -  End-shore  discharge  estimate;  ra’/s  or  ft'/s 

7  -  Ending  distance  (boat  to  shore);  w  or  ft 

8  -  Starting  depth  of  middle  layer  (or  ending  depth  of  top  layer);  m  or  ft 

9  -  Ending  depth  of  middle  layer  (or  starting  depth  of  bottom  layer);  m  or  ft 

1  NUMBER  OF  BINS  TO  FOLLOW 

2  MEASUREMENT  UNIT  -  cm  or  ft 

3  VELOCITY  REFERENCE  -  BT,  LAYER,  or  NONE 

4  INTENSITY  UNITS  -  dB  or  counts 

5  INTENSITY  SCALE  FACTOR  -  in  dB/count 

6  SOUND  ABSORPTION  FACTOR  -  in  dB/ra 


7 


1 

2 

3 

4 

5 

6 
7 

a 

9 

ID 

11 

12 

13 


DEPTH  -  Corresponds  to  depth  of  data  for  present  bin  (depth  cell); 
blanking  value;  in  m  or  ft. 

VELOCITY  MAGNITUDE 
VELOCITY  DIRECTION 

EAST  VELOCITY  COMPONENT  -  East(+)/Ucst(-) 

NORTH  VELOCITY  COH.pONENT  -  Horth(+)/Sauth(-) 

VERTICAL  VELOCITY  COMPONENT  -  Up(+)/Down(-) 

ERROR  VELOCITY 

ECHO  INTENSITY  -  Beam  1 

-  Beam  2 

-  Beam  3 

-  Beam  4 

PERCENT-GOOD 

DISCHARGE 


includes  ADCP  depth  and 


« 


APPENDIX  B 


ver  1.0  SOURCE  CODE 


1 


I 


c  a  title  subroutine 
call  clearscrccn 

writc(*,  *) '  Welcome  to  GPSQ  vci  1 .0  for  DOS' 

dok-1,  10 
write  C",  *) 
end  do 

write  (*,  *)  'press  F.NTI'R  to  continue' 
read  (*,  ’") 

call  tilcname 

end 


c- 


subroutine  filename 

c  a  subroutine  for  entering  adep  and  navigation  input  tile  names 
character‘*‘70  adep  file,  nav  file 
call  clearscreen 


writeC*,  *)  'ploase  enter  ADCP  ascii  file  name  and  path' 
write(*,  *) 

read(*,  '(a70)')  adep  file 

write(*, 

call  clearscreen 

writeC",  "')  'please  enter  NAVlGA'flON  file  name  and  path' 
write(*,  *) 

read(*,  '(a7())')  nav  file 
call  clearscreen 

call  data  out(adcp  file,  nav  file) 
end 


subroutine  data  out(adcp  file,  nav  file) 
c  a  subroutine  for  selecting  the  units  of  measurement  for  data  output 
real  units,  units2,  units3 
integer  option4,  error 


1 


cliaracter'''70  ft  cm*l,  adcp  file,  nav  file 
open  (4,  file  -  adcp  file,  status  =  'old') 
call  clearscreen 
1 0  continue 

writef*,  *)  'OPTIONS  FOR  DATA  OUTPU  T  .' 
writcC',  "') 

writcC*',  ’^)  '(1)  HNGLISH  units  (ft,  ll/s,  cfs)' 
writcC*',  *) 

write(*,  '")  '(2)  SI  units  (m,  m/s,  cuniccs)' 
doj=l,  10 
writeC',  *) 
enddo 

writeC",  *)  'please  enter  an  option  number' 
read  option4 

20  continue 

c  'units'  is  a  factor  used  to  convert  data  into  the  desired  outjiut 

select  case(option4) 
case( 1 ) 
units  3.281 
case(2) 
units  =  1 .0 
case  default 
call  clearscreen 

write  (*',  *)  'please  re-enter  an  option  number' 
goto  10 
end  select 


do  30,  j  1,  [] 
read  (4,  1 00) 

30  continue 
100  fomiat(70x) 

c  Adcp  data  file  is  read  to  determine  the  measurement  unit  'll  cm' 
c  'units2'  is  used  to  convert  the  ADCP  data  into  a  form  whereby 
c  'units'  can  be  applied 

c  'units3'  is  req'd  because  if  the  measurement  unit  is  SI  then  velocity 
c  data  is  given  in  cm  not  m 

read  (4,  '(3x,  al )',  end^40,  err -40,  iostat  error)  II  cm 

select  case(fl  cm) 

case('f) 

units2  =-  units/3.281 
units3  -  I 


I 


case('c') 

units2  "  units/ 1 .0 
units'!  Ofll 
end  select 

close  (4) 

call  q  mcthod(units,  units2,  units.l,  adcp  file,  nav  file) 

40  continue 
if(enoi  .nc.  0)then 
call  clearscreen 

write  'CANNOT  READ  measurement  unii' 
write  (*,♦)  'liRROR  in',  adcp  file 
stop 
endif 

end 


c  a  subroutine  foi  entering  preference  of  power  or  conslanl  method  for 
e  top  layer  discharge  estimation 

subroutine  qjnethod(units,  units2,  units3,  adcp  file,  nav  file) 

real  units,  units2,  units.1 

integer  option2 

character'* 70,  adcp  tile,  na^  file 
call  clearscreen 
10  continue 

write(*,  *)  'OPTIONS  FOR  DlSCHARGl-  CAI.CUEA'I  ION  ,' 
writc(*,  *) 

wrilcC*,  *)'(!)  use  "POWER"  method  for  top  layer  estimation' 
write(*,  '*) 

write!  ■*,  *)  '(2)  use  "CONS  I  AN  f"  method  for  top  layer  estimation' 
doj  1,  10 
writeC*,  ■*) 
enddo 

wi'ilcC*,  *)  'please  enter  an  option  number' 
read  ■*,  option2 

select  case  (option2) 
case( 1 ) 

call  discharge  1  (units,  units2,  units3,  adcp  iile,  nav  file) 


• 

case(2) 

call  discharge2(units,  units2,  iinits3,  adcp  iile,  nav  file) 
case  default 

call  cleai  scrccn  • 

write  (*,  *)  'please  re-enter  option  number' 
goto  10 
end  select 

end  » 


C' 


c  a  subroutine  for  calculating  discharge  using  the  [lowei  method 

subroutine  discharge  1  (units,  units2,  units3,  adcp  tile,  nav  file) 

real  dcpth_bin{100),  vcl  niag{10()),  vcl  dirn(l()0),  corr  w, 
i  c  vf.;i(]00),  n  vel(lOO), /.  vel(lOO),  db  avg(  100),  depth  avg, 

I  error  vel(10()),  total  binq,  toplaycr  q,  blnilayer  q, 

I  a  avg,  ensemble  q,  line  grad,  blank  depth,  a(  100), 

(  power  coeir,  bin  vel(lOO),  bin  q(  100),  adcp  range(  100), 

(  bin  size,  line  dim,  total  1,  total2,  total3,  total4,  bank  q, 
i  bin  velavg,  start  deptli,  start  vcl,  units,  units2,  units3, 

^  corr  e,  corr  n,  dist 

integer  ic,  ic2,  ic3,  num  ens,  ens  num,  id 

cliaracter''‘70  label  1,  label2,  adcp  file,  nav  file 

ic  ■  () 
ic2  I 
ic3  0 
tolal4  0 


call  velocity{units,  units2,  units3,  num  ens,  line  grad,  bin  size, 

(blank  depth,  adcp  file,  nav  iile) 

close  { I ) 

9 

open  ( 1 ,  file  'vclocity.dat',  status  'old') 
open  (2,  file  -  'gpscorrii.out',  status  'unknown') 

read(  1 ,  '(a70)')  label  1 

writc(2,  *)  label!  • 

read(l,'(a70)')labcl2 
wi  ite(2,  *)  labcl2 
read(l,  ’*') 


call  clcarscreen 

write(*,  *)  'PLEASK  ENTER  A  POWER  EAW  CO-EFFJCIENT' 

write  (*,  *) 

read  power  coeff 

call  clcarscreen 

write  (2,  *)  'top  layer  discharge  estimation  method:  power' 
write  (2,  *)  'power  co-efficient  power  coeff 
write  (2,  *)  'ensemble,  dist,  easting,  northing,  width,  depth,  bin 
^  q,  top  q,  btm  q,  total  q' 
write  (2,  '^)  'power  co-efficient power  coclT 

depth  bin(O)  ==  blank  depth  -  bin_size/2 

10  continue 
total  1  =  0 
total2  =  0 
total3  0 

20  continue 
ic  ic  t  1 

read  (1,  err-30,  end=30)  id,  ens  nuni,  dist,  corr  e, 

t  corr  n,  corr_w,  depth  avg,  depth  bin(ic),  vel_mag(ic), 

)  vel  d>nt(ic),  c  vel(ic),  n_vel(ic),  z  vel(ic),  error  vcl(ic), 

^  db_avg(ic),  adcp  range{ic) 
if  (id  ,eq.  ic2)  then 
goto  20 
endil' 
do  k^-'  1 ,  2 
backspace  ( 1 ) 
enddo 

read  (1,  err  30,  end=^30)  id,  ens  nuni,  dist,  corr  e,  corr_  n, 

I  ctnr  w,  depth  avg 

'a(  )'  is  equal  to  the  'non-adcp'  parameters  ir:  the  power  law  equation 

30  cr)ntinue 
ic  ic  -  I 

line  dim  ((atan(l/linc  giad)) 180/3.14159)  t  90 
do  40,  k= ! ,  ic 

bin  vel(k)  ^  vel  mag(k)’*‘cos((linc  dirn-vel  dim(k))*3. 14159/180) 
bin  q(k)  bin  vel(k)*corr  w*(depth  bin(k)  -  depth  bin(k-l)) 
a(k)  -  bin  vel(k)/( ''depth  avg  -  depth  bin(k))’'‘ ’“power  cocti) 
total  1  "=  bin  q(k)  t  total  1 
total2  =  a(k)  +  total2 
total3  bin  vel(k)  t  total3 
40  conlinuc 


total  bint]  total  1 
a  avg  total2/ic 
bin  velavg  totabVic 

toplayer  q  -  corr  w  *  blank  depth  *  a  avg  * 
t  ((depth  avg  -  blank  depth/2)* *powcr  coeft") 
bimlayer  q  -  a  avg  *  ((((depth  avg  - 
+adep  range(ic))**7)/2)**power  coefl')*corr  w 
ensemble  q  "  total  binc]  +  toplayer  q  t  btmlaycr  q 


write  (2,  1 00)  ens  num,  dist,  corr  e,  corr  n, 
( corr  w,  depth  avg,  total  binq,  toplayer  q, 
tbtmlayer  q,  ensemblc  q,  a  avg 
100  Ibrmat  (i5, 1^.2,  2fl3,2,  7110.2) 

if  (ic2  .cq.  1)  then 
start  depth  =  depth  avg 
start  vel  total3/ic 
endif 


total4  ensemble  q  ^  total4 
if  (ic.1  cq,  num  ens)  then 
goto  50 
endif 
ic  -  0 

ic2  -  ic2  1 
goto  10 
50  continue 

end  depth  ^  depth  avg 
end  vel  bin  velavg 

call  nearshore(bank  q,  stait  depth,  start  vel,  end  depth,  end  vel, 
^  units) 


call  clearscreen 
lolal4  ^  bank  q  h  total4 

if  (units  ,eq.  1 .0)  then 

write  (2,  *)  'total  discharge  ',  total4,  '(cumecs)' 
write(*,  *)  'TOTAL  DISCHARGE  - ',  total4,  '(cumecs)' 
do  k  1,  10 
write(*,  *) 
cnddo 
end  if 


it' (units  eq.  .1,28 1 )  then 
write  (2,  *)  'total  discharge  "" total4,  '(cfs)' 
writc(*,  *)  'TOTAL  DISCI lARGli  total4,  '(cfs)' 
do  k^l ,  10 
writeC",  *) 
enddo 
end  if 

60  continue 
end 


c  a  subroutine  for  calculating  discharge  using  the  constant  method 
c  a  power  law  fit  is  still  used  for  the  near-bed  discharge,  as  such 
e  'discharge2'  is  very  similar  to  'discharge  T 

subroutine  discharge2(unils,  units2,  units.l,  adcp  tile,  nav  tile) 

real  dcpth_bin(  1 00),  vel_  mag(100),  vel  dirn(UK)),  corr  w, 

-l  e  vel(lOO),  n  vcl(10()),  z  vcl(lOO),  db  avg(l()0),  depth  avg, 
(errc'r  vel(lOO),  total  binq,  toplayer  q,  btmlayer  q, 
t  a  avg  ensemble  q,  line  grad,  blank  depth,  a(IOO), 
t  power  coetf,  bm  vel(lOO),  bin  q(  100),  adcp  rangc(  100), 

I  bin  size,  line  dini,  total  1,  total2,  total!,  tolal4,  bank  q, 
t  bin  velavg,  vel  top,  start  deplh,  start  vel,  units,  uriits2, 

I  units.!,  corr  e,  corr  n,  dist 

integer  ic,  ic2,  id,  num  ens,  ens  num,  id 

character'''70  label  1,  label2,  adcp  tile,  nav  tile 

ie  0 
ic2  1 
id  0 
tolal4  0 

call  velocil,y(units,  units2,  units.!,  num  cn.s,  line  ^rad,  bin  size 
I  blank  depth,  adcp  tile,  nav  tile) 

close  ( 1 ) 

open  (1,  tile  'velocity.dat',  .status  'old') 
open  (2,  tile  'gpscorrq.out',  status  'unknown') 

read(  l,  '(a70)')  label  1 
write(2,  ♦)  label  1 
read(l,  '(a70)')  Iabel2 


sasts 


wriie(2,  *)  label2 
read( I,  *) 


depth_bin(0)  =  blank  depth  -  bin  size/2 
cal!  clearscreen 

write(*,  *)  'Please  enter  a  power  law  co-efficient ' 

write  C*,  ’") 

read  *,  power  coeff 

call  clearscreen 

write  (2,  *)  'top  layer  discharge  estimation  method:  constant' 
write  (2,  '*‘)  'power  coefficient  for  bottom  layer  =',  power  coetf 
write  (2,  *)  'ensemble,  dist,  easting,  northing,  width,  depth,  bin 
q,  top  q,  btm  q,  total  q' 


1 0  continue 
total  1  ""  0 
total2  =  0 
totals  0 

20  continue 
ic  ic  H  1 

read  (1,  en  -30,  cnd=30)  id,  ens  num,  dist,  corr  e,  corr  n. 

H  corr _w,  depth  avg,  depth  bin(ic),  vcl_mag(ic),  vcl_dim(ic), 

H  e  vel(ic),  n_vel(ic),  z_vel(ic),  crror_vcl(ic),  db  avg(ic), 

H  adcp__range(ic) 
if  (id  ,cq.  ic2)  then 
goto  20 
endif 
do  k  =  l,2 
backspace  ( 1 ) 
enddo 

read  (1,  *,  cn-30,  end-30)  id,  ens  num,  dist,  corr  e,  corr  n, 

I  corr  w,  depth  avg 

30  continue 
ic  ic  -  1 

line  dim  =  ((atan(l/line  jrad))*  180/3  14159)  H  90 

.1  ~  4  /■\  I  1  *.  .. 

UO  HU,  K  1,  le 

bin  vel(k)  vel  mag(k)''‘cos((line  dirn-vel_diin(k))*3, 141 59/1 80) 
bin  q(k)  ^  bin  vel(k)  ♦  corr  w  *  (depth  bin(k)-depth  l)in(k-l)) 
a(k)  "  bin  vel(k)/((depth  avg  -  depth  bin(k))'*‘*power  coeif) 
total  1  bin  q(k)  +  total  1 
total2  ""  a(k)  i  total2 
totals  -  bin  vel(k)  h  totals 

40  continue 


vel  top  -“=  vel  mag(l)*cos((line  dim  -  vel  dirn(l))*3. 141  59/180) 
ic3  ^  ic3  +  1 
totalbinq  =  total  1 
a  avg  =  total2/ic 
binvelavg  =  total3/ic 

to  player  q  =  corrw  *  blank  depth  *  vel  top 

btmlayer  q  =  a  avg  *  ((((depth  avg  - 

+adcp  range(ic))**7)/2)**power_coeff)  *  corr  w 

ensemble  q  =  total  binq  +  toplayer  q  +  btmlayer  q 

wnte  (2,  100)  ens  nutn,  dist,  corr  e,  corr  n,  corr  w,  depth  avg, 
+total  binq,  toplayer  q,  btmlayer  q,  ensemble  q,  a  avg 
100  format  (i5,  f8.2,  2fl3.2,  7fl0.2) 

if  (ic2  .eq.  1)  then 
start_depth  =  depthavg 
startvel  =  total3/ic 
endif 

total4  =  ensemble_q  +  total4 
if  (ic3  .eq.  num  ens)  then 
goto  50 
endif 
ic  =  0 

ic2  =  ic2  +  1 


goto  10 
50  continue 

end  depth  =  depth  avg 
end  vel  =  bin  velavg 

call  nearshore(bank  q,  start  depth,  start  vel,  end  deptli,  end_vel, 
+units) 

call  clearscreen 
total4  =  bank  q  +  total4 

if  (ui  4s  .eq.  1 .0)  then 

write  (2,  *)  'total  dischaige  =  total4,  '(cumecs)' 
w^te(^  *)  'TOTAL  DISCHARGE  =  total4,  '(cumecs)' 
do  k  -^1,  10 
write(’*‘,  *) 
enddo 
eii.  if 


if  (units  eq.  3.281)  then 
write  (2,  *)  'total  discharge  total4,  '(cfs)' 
write(*,  *)  'TOTAL  DISCHARGE  =  totalA  '(cfs)' 
dok=l,  10 
writeC*',  ’*') 
enddo 
end  if 

60  continue 
end 


c  a  subroutine  used  to  calculate  the  missing  discharge,  if  any,  at  the 
c  start  and  end  of  a  transect 

subroutine  nearshore(bank_q,  start  depth,  start  vei,  end  depth, 
+end  vel,  units) 

real  bank  q,  start  dist,  end_dist,  start  q,  end  q,  start  depth, 

H  start  vel,  end  depth,  end  vel,  units 

integer  options 

call  clearscreen 


05  continue 

write(*,  *)  'NEAR  SHORE  DISCHARGE  OPTIONS...' 
write(*, "') 

writeC^, "')  '(1)  CALCUi.ATE  near  shore  discharges' 
vvnte(*,  *) 

writeC",  *)  '(2)  DO  NOT  CALCULATE  near  shore  discharges' 
writeC*,  ’") 

writeC*, 'please  enter  an  option  number' 
read  options 

select  case(option3) 
case(  1 ) 
goto  10 
case(2) 
goto  20 
case  default 
call  clearscreen 

write  C',  *)  'please  re-enter  option  number' 
goto  05 
end  select 


1 0  continue 


call  clearscrcen 


if  (units  ,eq.  1.0)  then 

write(*,  *)  'enter  distance  to  bank  at  S'I'AR’l'  of  transect  (m)' 
read  *,  start  dist 
writeC**,  *) 

write(*,  *)  'enter  distance  to  bank  at  HND  of  transect  (m)' 
read  *,  end  dist 
end  if 

if  (units  ,eq.  3.281)  then 

\vrite(*,  *)  'enter  distance  to  bank  at  START  of  transect  (ft)' 
read  *,  start  dist 
writeC",  *) 

writeC",  *)  'enter  distance  to  bank  at  END  of  transect  (ft)' 
read  end  dist 
end  if 

call  clearscrcen 

start  q  0.707  start  depth  start  dist  ■"  start  vcl/2 

end  q  -  0.707  *  end  depth  *  end  dist  end  vel/2 

bank  q  start  q  +  end  q 

goto  30 

20  continue 
bank  q  =  0 

30  continue 

end 

c - - - 

subroutine  velocity(units,  unit.s2,  units3,  num  ens,  line  grad, 

H  bin  size,  blank  depth,  adcp  file,  nav  file) 

c  a  subroutine  for  combining  the  vessel  velocity  and  the  adcp  measured 
c  water  velocities 

real  depth(  1 00),  mag  adept  1 00),  dim  adcp(  1 00),  e  adcp(  1 00), 

I  n  adcp(  1 00),  z  adcp(  1 00),  error  vel(  1 00),  db  1  ( 1 00),  db  2(  1 00), 
+db  3(100),  db_4(  100),  per  j50od(  100),  e  boat,  n  boat,  corr  c, 

I  corr  n,  corr  w,  depth  nvg,  dist,  c_water(  100),  n  waler(  100), 

I  z  water(  1 00),  mag  ■water(  1 00),  dim  watcr(  1 00),  db  avg(  1 00), 
t  junk  1  ( 1 00),  junk2(  1 00),  junk3(  1 00),  line  grad,  bin  size, 


+adcp  angle,  adcp  range,  bin  ratio,  blank  depth,  mag  dccln, 

^  blank  dist,  units,  units2,  units3,  jutik4 

integer  ens  num,  good  bins,  nuinbins,  enseniblenum,  num  ens,  ic, 
<  ic2,  error,  error2,  error3,  option  1 

character*70  label  1,  label2,  adcp  file,  nav  file 

ic  =  0 
ic2  -0 

call  clearscrccn 
05  continue 

write(*,  ♦)  'DEPTH  CALCULATION  OPTIONS  ..  ' 
writcC",  "') 

writcC**,  '*')  '(1)  use  ADCP  beam  depths' 
writeC*',  *) 

writeC^,  *)  '(2)  use  FATHOMLTI-R  data' 
doj- 1,  10 
write(*,  ’*') 
enddo 

’writc(*,  'please  enter  an  option  number' 
read  option  1 

select  case  (option  1 ) 
case  (1) 

call  ens  headerl  (units,  num  ens,  line  grad,  nav  file) 
case  (2) 

call  cns_headcr2(units,  num  ens,  line  grad,  nav  file) 
case  default 
call  clearscreen 

write  (*,  *)  'please  re-enter  option  number' 
goto  05 
end  select 

rewiiul  (3) 

open  (1,  file  'vclocity.dat',  status  --  'unknown') 
open  (4,  file  =  adcp  file,  status  =  'old') 

cal!  clearscreen 

write(*,  ''‘)  'please  enter  ADCIP  HHAM  ANGLE  (deg)' 

wi  ite  (*,  *) 

read  *,  adcp  angle 

call  clearscreen 

write(*,  '*')'irNOT  already  collected.,.' 

write(*,  *)  'please  enter  the  local  magnetic  declination  (deg)' 

write  (*,  *) 


read  *,  mag  decln 

call  clearscreen 

read(4,  '(a70)')  label  1 
read(4,  '(a70)')  label2 
write(l,  *)  label  1 
write(l,  *)  label2 

write(  l,  100)  'id,  ens,  dist,  east,  nortli,  width,  depth,  bin,  vcl 
“i  mag,  vel  dim,  east  vel,  north  vel,  vert  vel,  en  or  vel,  db  avg,  a 
+dcp_range' 

100  format(al20) 

read(4,  *,  err=80,  end^95,  iostat=enor)  bin  size,  junk4, 
t  adcpdepth 

bin  size  -  (bin  size/ 100)  *  units 
adcp  depth  ^  (adcp  depth/ 100)  units 

10  continue 
ic  “  0 

20  continue 
ic2  =  ic2  H-  1 

read  (4,  eri~80,  end=95,  iostat-error)  (junkl(i),  i-^l,  7), 

-t-ensemblenum,  (juriik2(i),  i=l,  13),  depth  bl,  depth  b2,  depth  b3, 
+dcpth_b4,  (junk3(i),  i  1,  19),  nurnbins 
depth  avg  =  (dep'.h  bl  -t  depth  b2  +  depth  b3  +  depth  b4)/4 
depth  avg  --  dopth  avg  ♦  units2 

30  continue 
if  (option  1  .eq.  1)  then 

read  (3,  *,  err =85,  end  95,  iostat  error!)  ens  mini,  dist, 

I  corr  e,  corr  n,  corr  w,  e  boat,  n  boat 
cndif 

if  (optionl  .eq.  2)  then 

read  (3,  err--85,  end=95,  iostat  error!)  ens  nuin,  dist, 

^  corr  e,  corr  n,  corr  w,  e  boat,  n  boat,  depth  avg 

if  (ens  nuni  .It.  enseinblenum)  then 
goto  30 
endif 


40  continue 
ic  ic  )  1 

read  (4,  *,  end=50,  err=90,  iostat=error3)  dcpth(ic), 

I  mag  adcp(ic),  dim  adcp(ic),  c  adcp(ic),  n  adcp(ic),  z  adcp(ic). 


terror  vel(ic),  db  l(ic),  db  2(ic),  db  .^(ic),  db  d(ic), 

t  per  ^ood(ic) 

ir(ic  ge.  numbins)  goto  50 

goto  40 

50  continue 

if(cns  num  .gt.  ensemblenum)  then 
goto  20 
cndif 

do  j  1,  ic 

dcpth(j)  dcptli(j)  *  unils2 
niag  adcpO)  mag  adcp(j)  *  units2  *  units3 
dim  adcp(j)  -  dim  adcp(j)  -  mag  decln 
il^dirn  adcp(j)  gt,  300)  ilten 
dim  adep(j)  ^  dim  adcp(j)  -  360 
endil’ 

/,  adcp(j)  ^  /  adcp(j)  *  iinits2  *  units3 
error  vcl(j)  --  error  vel(j)  *  units2  *  units3 
cnddo 

do  60,  k  1 ,  ic 

c  adcp(k)  ^  inag  adcp(k)  *  siri((dirn  adcp(k))*3. 14159/180) 
n  adcp(k)  mag  adcp(k)  *  co!i((dirn  adcp(k))*3. 14159/180) 

0  watcr(k)  ^  c  boat  ^  c  adcp(k) 
n  walei(k)  ii  boat  i  ii  adcp(k) 
z  water(k)  -  z  adop(k) 

mag  watcr(k)  "  (c  walcr(k)'''*2  i  n  walor(k)*'''2)**0.5 
dim  watcr(k)  (aoos{n  water(k)/inag  watcr(k)))*  180/3, 14159 
db  avg(k)  (db  l(k)^db  2(k)^db  3(k)idb  4(k))/4 
60  continue 


adep  range  (depth  avg  -  adep  deptli)*(eos((adep  angle)* 

1,3,14159/180))  I  adep  depth 

blank  dist  depth(l)-bin  size/2  -  adej)  depth 

blank  depth  -  adep  depth  l  blank  (list 

bin  range  adep  range  -  blank  depth 

bin  ratio  bin  rangc/l)in  size 

good  bins  aint{bin  ratio) 

do  70,  k  1 ,  good  bins 
ir(mag  water(k)  gt,  95)  then 
goto  70 
endil' 

write  ( 1 ,  200)  ic2,  ens  nuin,  dist,  eon  e,  con  n,  eon  w, 

I  deptli  avg,  depth(k),  mag  water(k),  dim  water(k),  e  watei{k), 
t  n  water(k),  z  water(k),  error  vel(k),  db  avg(k),  tidcp  lange 
70  continue 


H 


200  format  (2i5,  f8,2,  2n3,2,  1  ino.2) 
goto  10 

80  continue 
if  (error  ,ne,  0)  then 

write(*,  *)  'error  in  adcp  header,  ensemble  no,',  ic2 
stop 
end  if 


85  continue 
if  (error2  .  ne.  0)  then 

write(*,  '*)  'error  in  navfile,  ensemble  no.',  ic2 
stop 
endif 

90  continue 
if  (error3  gt.  0)  then 

writeC",  'error  in  adcp  bin  data,  ensemble  no.',  ic2 
stop 
endif 


95  continue 
end 


c- 


subroutine  ens  header  1  (units,  num  ens,  linc  jjrad,  nav  j'ile) 

c  a  subroutine  for  calculating  the  survey  line  definintion,  corrected 
c  ensemble  positions  and  widths,  and  vessel  velocities 

real  e  end(lOO),  n  cnd(lOO),  e  vel,  n  vcl,  c  start,  n  start, 
t  line  grad,  linejnt,  corr  e,  coit  n,  corr  w,  sum  e, 

I  sum  esq,  sum  n,  sum  en,  con  estait,  corr  nstart,  last  c,  last  n, 
t  corr  eavg,  corr  navg,  units,  dist 

integer  ic,  ens  num(10()),  num  ens,  time(lOO),  time  start 

character''‘70  nav  file 

sum  e  "  0 
sum  esq  =  0 
sum  n  "  0 
sum  en  “  0 
last  e  0 
last  n  0 
ic  “  0 


call  read  gpsl(num  ens,  time  start,  e  start  n  start,  time, 

+eiis  num,  e  end,  n  end,  nav  file) 

open  (3,  status  =  'scratch') 

suin_c  =  e  start  +  sum_e 

sum  esq  ^  e_start’'"'‘2  -i  sum  esq 

sum  n  "  n  start  t-  sum  n 

sum  en  =  e  start’''n  start  +  sum  en 

ic  =  num_ens  +  ic 
do  10j=l,  ic 

sume  =  sum_e  +  c_cnd(j) 
sumesq  =  sum_esq  +  e_end(j)*'*'2 
sum  n  =  sum  n  +  n  end(j) 
sum_en  =  sumen  +  e_end(j)'*'nerid(j) 

1 0  continue 

ic  ~  ic  t  1 

line^rad  =  (ic'''sum  en  -  sum  j;*suni  n)/(ic*sum  esq  -  sum  e*''‘2) 
linejnt  -  (sum_esq*surn  n  -  sum_e’‘'sum_en)/(ic'*sum_csq  -  sum_e*'''2) 
ic  =  ic  -  1 

e  vel  -  (c_  end(l)  -  e_start)'*'units/(time(  1)  -  time  stait) 
n_vel  =  (n_end(l )  -  n  start) '''units/(timc(  I )  -  time  start) 

corr  estart  =  (line_grad’'‘n  start  -i  e  start  -  iinej»rad*lincjnt) 

-I /(line  grad* '''2  +  1) 

corr  nstart  “  line  grad*corr  estart  +  line  int 

corr  c  (line  grad*n_end(1)  t  c  cnd(l)  -  line  grad*line  int) 
+/(linc_grad**2  +  1) 
corr  n  =  line  grad  *  corr  e  i  lincint 

corr  w  (((corr  e  -  con  estart)**2  l  (corr  n  -  corr  iistait)**2) 
+**0.5)*units 

dist  ^  corr_w/2 

corr  eavg  =  (corr  estail  +  corr  e)/2 
corr  navg  =  (corr  nstart  t  corr  n)/2 

write(3,  '*')  ens  num(l),  dist,  corr  eavg,  corr  navg,  corr  w,  e  vel, 

+n  vel 

last  e  corr  e 
last  n  =  con  n 


do  20,  j-2,  ic 


e  vel  ^  (e  cnd(j)  -  c  cnd(j-l))''  units/(tiiiic(j)  -  limc(j-l)) 
n  vcl  (n  cnd(J)  -  n  end(J-l))*iinits/(timc(J)  -  timcO-l )) 

corr  e  =  (linc^rad*n  end(j) -I  e  cnd(j)  -  line  {•rad*linc  int) 
^  /(line  grad**2  1) 
coir  n  line  grad*corr  e  )  line  int 


coir  w  "  (((corr  e  -  last  e)*''2  h  (corr  n  -  last  n)**2)**().5)* 
^  units 

dist  (((corr  e  -  corr  estart)*'*'2  l  (corr  n- corr  nstart)**2) 

^  **0.5)*units '  corr  w/2 


corr  eavg  (corr  e  ^  last  c)/2 
corr  navg  =  (con  n  t  last  n)/2 

write  (3,  *)  ens  tiuni(j),  dist,  corr  eavg,  corr  nuvg,  corr  w, 
‘  c  vel,  n  yel 

last  e  -  corr  c 
last  n  ■  corr  n 
20  continue 

end 


c- 


subroutine  ens  header2(units,  nuiii  ens,  line  grad,  nav  file) 

c  a  subroutii.e  lor  calculating  the  suiwey  line  deliniiition,  corrected 
c  ensemble  positions,  widths,  depths  and  vessel  velocities 


reale  eiid(10()),  n  end(lOO),  e  vel,  n  vcl,  e  stait,  ii  start, 

I  line  grad,  line  int,  con  e,  corr  ii,  corr  w,  sum  c, 

I  sum  esq,  sum  n,  sum  en,  corr  estart,  con  nstart,  last  e,  last  n, 
I  corr  eavg,  corr  navg,  units,  dist,  depth  avgdOO) 

integer  ic,  ens  num(  100),  nuni  ens,  timc(IOO),  time  start 

charactei  ’*'70  nav  tile 

sum  c  0 
sum  CSC]  0 
sum  n  0 
sum  cn  -  0 
last  c  () 
last  n  0 
ic  -  0 


call  read  gps2(num  cns,  time  start,  c  start,  u  start,  time, 
l  ens  nuin,  e  end,  n  end,  depth  avg,  nav  file) 

open  (.1,  status  'scratch') 

sum  c  e  start  i  sum  e 

sum  esq  c  start’"  *2  i  sum  esq 

sum  n  n  start  i  sum  n 

sum  en  e  start ’"n  start  i  sum  en 

ic  "  num  cns  h  ic 
do  1 0  j  1 ,  ic 

sum  e  sum  c  l  e  end(j) 

sum  esq  sum  esq  i  e  end(j)*’*‘2 

sum  n  sum  n  h  ii  end(j) 

sum  CM  sum  en  i  e  end(j)'''n  end(j) 

10  continue 

ii:  ic  I  1 

line  grad  (ic^sum  en  -  sum  e’''sum  n)/(ic'"sum  esq  -  sum  e''""2) 
line  int  (sum  esq*sum  n  -  sum  e'"sum  en)/(ic*sum  esq  -  sum  e* *2) 
ie  ic  ■■  1 

do  15,  j  1,  ic 

depth  avg(l)  depth  avg(j)  *  units 
15  continue 


e  vel  (e  end(l)-e  starl)*uiiits/(time(  1 )  -  lime  start) 
n  vel  (n  end(l)-n  start )'"unit!;/(time(  I )- time  start) 

ernr  eslart  tline  giad""!!  start  l  e  start  line  grad*Iine  int) 

I /(line  grad'"  "2  i  1) 

con  nstari  line  grad^corr  eslart  i  line  int 

corr  e  (line  giad*n  cnd(l)  I  e  end(l)-liiic  grad'''line  ml) 

I /(line  grad’"'''2  i  1) 

corr  n  line  grad*corr  e  )  line  int 

corr  w  (((con  e  -  corr  eslarl)’'‘*2  i  (corr  n  -  corr  nstarl)*’"2) 
I  **()  5)*unils 

(list  corr  w/2 

corr  eavg  (corr  estarl  (  corr  e)/2 
corr  navg  (con  nstart  i  corr  n)/2 


write(5,  *)  ens  num(  1 ),  dist,  corr  eavg,  con  iiavg,  con  w,  c  vel. 


I n  vel,  depth.  avg(  1 ) 

last  e  corr  e 
last  n  corr  n 

do  20,  j  2,  ic 

e  vel  (e  end(j)  -  e  end(j-l  ))*units/(timc(j)  -  limc(j-l)) 
n  vel  (n  end(J)  -  n  end(J-l))*units/(time(J)  -  time(j-l)) 

corr  e  -  (line  ^rad*n  end(j)  t  e  cnd(j)  -  line  grad*line  int) 
I  /(line  grad**2  i  1) 
corr  n  ^  line  j5rad*corr  c  ■*  line  int 


corr  w  -  (((corr  c  -  last  c)''”''2  +  (corr  n  -  last  n)*’''2)**0.5)* 
t  units 

disf  (((corr  c  -  corr_estart)**2  h  (corr  n  -  corr  nstai1)**2) 

I  '*'*().5)'''units  -  corr  w/2 

corr  eavg  "  (corr  c  I  last  e)/2 
corr  navg  ■  (corr  n  H  last  n)/2 

write  (3,  *)  cns_runi(j),  dist,  corr  eavg,  corr  iiavg,  corr  w, 
t-  c  vcl,  n  vel,  de[jth  avg(j) 

last  c  corr  e 
last  n  -  con  n 
20  continue 

end 


subroutine  read  gpsl(ic4,  time  start,  e  start,  n  start,  • 

t  ens  time,  ens  luim,  e  end,  n  end,  nav  tile) 

c  a  subroutine  for  extracting  ensemble  data  from  gps  file  which 
c  does  not  include  depth  data 

a 

reale  Start,  it  start,  e  cnd(lOO),  n  ciid(lOO),  gps  e(lOOO), 

^gps  n(10(J0) 

integer  time  start,  hour  .start,  min  start,  sec  start, 

l  ens  time(  100),  gps  hour(  1000),  gps  rnin(  1000),  gps  scc(  1000),  ic,  • 

I  ic4,  error,  ens  num(lOO) 


character*  1  chk,  nav  nie*70 


ic  =  0 
ic4  •=  0 

open  (5,  tile  =  nav  file,  status  =^- old') 

01  continue 
read  (5,  '(al)')  chk 
select  case(chk) 
case('E') 
goto  02 
case('e’) 
goto  02 
case  default 
goto  0 1 
end  select 

02  continue 

read  (5,  200,  en  ■20,  end  20,  iostat-  eiror)  hour  start,  min  stail, 
+sec_stait,  e  start,  n  start 

titne_start  =  hour  start "'3 600  t  min_start'*'6()  t  sec  start 


do  j  1,2 
backspace  5 

enddo  • 

05  continue 
ic4  ■  ic4  t  1 

read  (5,  100,  crr--^40,  end'  40)  ens  nuin(ic4) 

10  continue  • 

ic  ic  H  1 

read  (5,  200,  err- -30,  end  -  30,  iostat  -erroi )  gps  hour(ic), 
tgps  min(ic),  gps  scc(ic),  gps  e(ic),  gps  n(ic) 

goto  10  • 

100  foimat(10x,  i4) 

200  format(i2,  lx,  i2,  lx,  i2,  fl2,2,  fl0,2) 

20  continue  • 

if  (error  ne.  0)  tlien 
writeC*",  *)  'error  in  gps  data' 
stop 
endif 

30  continue 
if  (error  ne.  0)  then 
ic  ic  2 

ens  tiiiie(ic4)  gps  hour(ic)'''3600  t  gps  rnin(lc)*60  I  gps  seo(ic) 


iHcns  tiine(ic4)  cq.  0)  then 

e  end(ic4)  ^  c  cnd(ic4-l)  i  (e  end(ic4-l)  -  c  end(ic4  2)) 
n  cnd(ic4)  n  cnd(ic4-l)  ^  (n  end(ic4-l)-n  end{ic4-2)) 
ons  tiine(ic4)  -  2*ens  time(ic4-l)  -  ens  timc(ic4-2) 
wntc(*,  *)  'error  at  ens  nutii',  ens  num(ic4) 
rcad(*,  *) 
ic4  ic4  )  1 
goto  40 
endif 

e  end(ie4)  gps  c(ic) 
n  end(ie4)  ^  gps  n(ie) 
baek space  5 
goto  05 
endif 

40  continue 
ic4  ic4  -  1 
end 


subroutine  read  gps2(ic4,  time  start,  e  stall,  n  start, 

H  ens  time,  ens  nuni,  e  end,  n  end,  deptli  avg,  nav  Itle) 

c  a  subroutine  for  extracting  ensemble  data  from  gps  lile  whieli  includes 
c  depth  data 

reale  start,  n  start,  e  end(  100),  n  end{  100),  gjis  e(l<)b()), 
igps  n(  1000),  depth  total,  depth  avg(  100),  dcpthflOOO) 

integer  time  start,  hour  start,  min  start,  sec  start, 
t  ens  time(  1 00),  gps  hour(  1 000),  gps  min(  1 000),  gps  see(  1 000),  ie, 
I  ie2,  icl,  ie4,  error,  ens  num(  100) 

character''' 1  chk,  nav  1ile'''70 

ic  0 
ic2  -  0 
ici  0 
ic4  0 

open  (5,  lile  nav  lile,  status  'old') 

1)1  continue 
read  (5,  '(a  I )')  chk 
selecl  case(cbk) 
ease('lf) 


goto  02 
case('e') 
goto  02 
case  default 
goto  0 1 
end  select 

02  continue 
ic2  ic2  1“  1 

read  (5,  200.  err"=20,  end=20,  iostat==error)  hour  start,  min  start, 
+SCC  start,  n  start,  c  start 

time  start  =  hour_start*3600  t  min_start*60  i  sec  start 

doj=l,2 
backspace  5 
enddo 

05  continue 
ic4  ~  ic4  I  1 

read  (5,  100,  err -^40,  cnd"40)  ens  num(ic4) 

10  continue 
ic  ■■  ic  t  1 

read  (5,  200,  err-30,  end-30,  iostat-=error)  gps  hour(ic), 
t  gps  niin(ic),  gps  sec(ic),  gps  _n(ic),  gps  e(ic),  depth(ic) 

goto  10 

100  format(  lOx.  i4) 

200  rormat(i2,  lx,  i2,  lx,  i2,  02,2,  fl0,2,  00,2) 

20  continue 
if  (error  ,ne,  0)  then 
write(*,  *)  'erroi'  in  gps  data' 
stop 
endif 

30  continue 
if  (error  ,ne,  0)  then 
ic  =  ic  -  2 
depth  total  0 

do  j  =  ic2,  ic 

depth  total  -  depth(j)  +  depth  total 
ic3  ic3  i  1 
enddo 

depth  avg(ic4)  ^  depth  total/ic3 

ens  time(ic4)  gps  hour(ic)'*' 3600  i  gps  min(ic)’'‘60  ^  gps  scc(ic) 


e  end(ic4)  gps  e(ic) 
n  end(ic4)  -  gps  n(ic) 


ic2  -  ic  +  1 
ic3  -  0 
backspace  5 
goto  05 
endif 

40  continue 


ic4  =  ic4  -  1 


end 


subrouiine  clearscreen 

do  k"l,  25 
write(*,  *) 
enddo 

end 


THE.END..»—{ 


c- 

c 

C' 


APPENDIX  C 


ADCPFILL  SOURCE  CODE 


c  a  subroutine  for  filling  in  bad  ensemble  data 

real  bdepthl(lO).  bdepth2(l()),  bdepth3(10),  heada  1(20), 

+heada  2(20),  heada  3(20),  headb  1(30),  headb  2(30),  headb  3(30), 
Hopbtm  1  ( 1 0),  topbtm2(10),  topbtm3(10) 

real  deptb](100),  depth2(100),  depth3(100),  veljnag  1  ( 1 00), 

+vel  mag2(100),  vel  jnag3(100),  vel  dirnl(lOO),  vcl  dirn2(100), 

•  vcl  dim3(  1 00),  east  vel  1  ( 1 00),  east  vel2(  1 00),  east_vel3(  1 00), 
+north_vel  1(100),  north  vel2(  100),  north  vel3(IO0),  veil  vel  1(1 00), 
3  vcit  vel2(  1 00),  vert_vel3(  1 00),  error  vel  1  ( 1 00),  error  vel2(  1 00), 
+error_vcl3(100),  db  j  1(100),  db  12(100),  db  13(100)^  db  21(100), 
3  db  22(100),  last  binq,  q  dill', 

3-db”23(100),  db_31(100),  db  32(100),  db  33(100),  db  41(100), 

+db  42(100),  db  43(100),  disehargcl(  100),  discharge2(  1 00), 
+dischargc3(100),  last  t,  dt,  bin  size,  units,  e  total, 

3  n  total,  bin  total,  1  total,  dc,  dn,  dl,  dbin  q,  dtop  q,  dbtm  q 

integer  time  1(10),  timc2(  10),  time3(  10),  ens  numl,  ens  num2, 

3-ens  num3,  ic,  ic2,  ic3,  ie4,  num  bins,  per  goodl(lOO), 

3  per  good2(  100),  per  good3(10()) 

character*70  label!,  Iabel2,  labcl3,  labcl4*50,  lllc  name,  11  cm*l 

ic  -  0 
ic2  =  0 
last  binq  ^  1 
last  top  0 
last  btm  --  0 
de  0 
dn  0 
dl  0 
dbin  q  0 
dtop  q  0 
dbtm  q  0 

call  clearscreen 

writc(*,  *) '  Welcome  to  AIX’I'FILI,  for  DOS' 

do  k-  1,  10 
writeC",  *) 
enddo 

write(*,  *)  'press  liN'fbR  to  continue' 
rcad(*,*) 

05  continue 
call  clearscreen 

wiitet*,  *)  'please  enter  adcp  ascii  lilo  name  and  i)alb' 
rcad(*,  '(a50)',  end  05,  err  05)  lile  name 


call  clcarscreen 

write(*,  *)  'please  wait  while  ADCPh'ILl.  checks...' 

write(*,  *) '  file  name 

do  ! ,  10 

write{*,  *) 

enddo 

open(  1 ,  file  ^  file  name,  status  =  'old') 
open(2,  file  =  'clean,  out',  status  =  'unknown') 
open(3,  status  =  'scratch') 

read(l,  '(a70)',  end=90,  eiT-==90)  label  I 
write(2,  ’•')  label  1 

read(l,  '(a70)',  end-90,  err=90)  label2 
wiite(2,  '•')  label2 
read(l, '")  bin  size 
backspace  1 

read(l,  '(a70)',  end=90,  err=90)  labcl3 
write(2,  *)  labels 


1 0  continue 

read(l,  end  =-90,  err=90)  (timcl(i),  i^l,  7), 
t  ens  numl,  (heada_l(i),  i=l,  13),  (bdepthl(i),  i=^l,  4), 
+(headb_l(i),  i=l,  19),  num  bins 


backspace  1 

read(  1 , 100,  end=90,  err=90)  ft  cm 
100  format(3x,  al) 
select  case(tt_cm) 
case('f) 
units  =  3.28 1 
case('c') 
units  ^  1 .0 
end:;clect 
backspace  1 

read(  1 ,  '(a50)',  end^"90,  err  90)  label4 


if  (bdepth  1(1)  .It.  i  .0)  then 

bdcpthl(l)  -  (bdepth  1(2)  +  bdepth  1(3)  +  bdepth l(4))/3 
endif 

if  (bdepth  1(2)  .It.  1.0)  then 

bdepthl(2)  =  (bdepthl(l)  i  bdepthl(3)  ^  bdepthl(4))/3 
endif 

if(bdcpthl(3)  It.  1.0)  then 

bdcpthl(3)  ~  (bdepthl(])  +  bdcpthl(2)  ^  bdcpthl(4))/3 
endif 

if  (bdepth  1(4)  ,  It.  1.0)  then 


1 


bdepthl(4)  =  (bdcpthl^i)  +  bdepthl(2)  )  bdcpthl(.'^))/3 
endif 


last  t  =  headb  1(2) 

20  continue 
ic  ==  ic  1 

read(l,  *,  cnd==00,  err^90)  depthl(ic),  vel  magl(ic), 
+vel_dirnl(ic),  east_vell(ic),  north_vell(ic),  vert_vell(ic), 
+error_vell(ic),  db  1  l(ic),  db  21(ic),  db  31(ic),  db  41(ic), 
+per  goodl(ic),  discharge l(ic) 
it(ic  .ge.  num  bins)  goto  30 
goto  20 

30  continue 

q  diff-  lieadb  1(1 1)  -  last  binq 
it^q  dift'.eq.  0)  then 
ic  =  0 
goto  10 
endif 


ic  =  ic  -  1 
35  continue 

write(2,  200)  (timel(i),  i=^l,  7),  ens  numl,  (heada_l(i),  i=l,  13), 
+(bdepthl(i),  i=l,  4),  (headb  l(i),  i=l,  19),  label4 
do  k=l,  num_bins 

write(2,  300)depthl(k),  vel  magi (k),  vcl  dirnl(k),  east  yell (k), 
+  north  vel  1  (k),  vert  vel  1  (k),  eit  or_vcl  1  (k),  db  l  1  (k),  db  2 1  (k), 
+  db  31(k),  db  41(k),  per_good  1  (k),  discharge l(k) 
enddo 

200  format(7i4,  i6,  518.2,  /,  8110  2,  417.2,  /,  511 0.2,  /,  5fl0.2,  /, 

+  3n4.2,  618.2,  /,  a50) 

300  fbrmat(r9, 1,  n  )  .2,  fl  1.3,  819. 1,  i7, 133.2) 
rewind  (3) 

40  continue 

last  binq  headb_l(l  1) 

topbtml(l)  "  headb  1(12)  -  last  top 

topbtml(2)  =  headb_l(13)  -  last  btm 

last  top  =  headb  1(12) 

last  btm  =  headb_l(13) 

read(l,  *,  end“90,  err--90)  (time2(i),  i^l,  7), 

+ens  nurn2,  (heada  2(i),  i^  l,  13),  (bdepth2(i),  i=l,  4), 

^  (headb_2(i),  i'^1,  19),  num  bins 

if  (bdepth2(  1 )  It.  1 .0)  then 

bdepth2(l)  =■  (bdcpth2(2)  H  bdepth2(3)  (  bdepth2(4))/3 


BMliililillinWiftlifri 


cndil' 

ii'(b(.lcptli2(2)  It.  1 .0)  then 

bdepth2(2)  (bdcpth2(  I )  i  bdcptli2(3)  i  l)deptb2(4))/3 
endir 

ir(bdeptb2(4)  It,  '.())then 

bdepth2(3)  (bdepth2(l )  )  bdcptli2(2)  l  bdeptli2(4))/3 
cndil’ 

ir(bdcptl]2(4)  ,lt.  1.0)  then 

bdeptli2(4)  (bdepth2(l)  i  bdcpth2(2)  i  bdepth2(3))/3 
endii' 

q  dill'  headb  2(11) -last  binq  i  dbin  q 
if(q  diir.nc.  0)  tlicn 
goto  00 
endir 

ie2  ie2  i  1 

dt  headb  2(2)  -  last  t 
last  t  headb  2(2) 

writc(3,  *)  (liinc2(i),  i  1,7),  headb  2(2),  dt 

do  50,  k  l,  mini  bins 
rcad(l,  400,  end  00,  err  00) 

50  eontinuo 
400  rornial(70x) 

goto  40 

00  eontinuc 

ir(ie2  gl.  0)  goto  70 
lime  I  time2 
ens  iiiiml  ens  imni2 
heada  I  heada  2 
bdepthl  bdeplh2 
headb  I  headb  2 
last  1  headb  1(2) 

headb  1(1)  headb  1(1)  l  dl 

headb  1(3)  headb  1(3)  i  dn 

headb  1(4)  headb  1(4)  i  de 

headb  1(5)  (head!)  I(3)**2  i  headb  1(4)**2)**0.5 

headb  1(11)  headb  1(11)  I  dbin  (| 

headb  1(12)  headb  1(12)  i  dtop  (| 

headb  1(13)  headb  1(13)  i  dbim  ii 


ie  0 


70  continue 
call  clearscreen 
ic  ""  0 

bin  total  =  headb  1(11) 
n  total  hcadb_l(3) 
e  total  headb  1(4) 

1  total  "headb  1(1) 

topbtm2(  1 )  =  headb_2(  1 2)  -  headb  1(12) 
topbtm2(2)  headb  2(13)  -  headb  1(13) 
rewind  (3) 

80  continue 
ic  =  ic  I  1 

read(l,  *,  cnd-90,  en-90)  depth2(ic),  vel  niag2(ic), 
t  vel  diin2(ic),  cas^  _vel2(ic),  north  vel2(ic),  vci1_vcl2(ic), 
t  error  vcl2(ic),  db_  I2(ie),  db_22(ic),  db  32(ic),  db  42(ic), 
t  per  good2(ic),  di  schargc2(  ic) 
it(ic  .It.  num  bins)  then 
goto  80 
endii’ 

ic3  =-■=  io2 
ic4  I 

do  85  k  I,  ic2 

topbtni3  (topbtml  *ic3  H- topbtin2*io4)/(ic3t  ic4) 
ens  nuni3  (ens  numl  *ic3  i  ens  nuin2*ic4)/(ic33-ic4) 
heada  3  (hcada  l*ic3  t  licada  2'*ic4)/(ic3 1  ic4) 
bdepth  -  (bdeplhl  *ic3  (  bdepth2*ic4)/(ic3  t  ic4) 
headb  <  (headb  I  '*'ic3  t  headb  2*ic4)/(ic3 1  ic4) 
depth3  depth  I 

vel  niag3  ^  (vel  magi *ic3  l  vd  mag2’''ic4)/(ic3  i  ic4) 

vel  dirn3  ■  (vcl  dirnl*ic3  i  vel  dirn2*ic4)/(ie3+ic4) 

ea.st  vel3  (east  vell*ic3  i  east  vel2*ic4)/(ic3  t  ic4) 

north  vel3  (north  vel  I  *ic3  t  north  vcl2*ic4)/(ic3  t  ic4) 

vert  vel3  (vert  vell*ic3  +  vert_vcl2’''ic4)/(ic3  I  ic4) 

error  vel3  (error  vell*ie3  l  error  vel2*ic4)/(ic3^-ie4) 

db  I3  (db  ll^ic^  t  db  I2*ic4)/(ic3  l  ic4) 

db  23  (db  2l*ie3  i  db  22*ic4)/(ic3-tic4) 

db  33  (db  3l*ic3  t  db  32*ic4)/(ic3  i  ic4) 

db  43  -  (db  4l*!c3  H  db  42*ic4)/(ic3t  ic4) 

per  good3  -  (per  goodl '''ic3  I  per  good2*ic4)/(ic3+ic4) 

re'id(3,  *)  (tinte3(ik),  ik  d,  7),  headb  3(2),  dt 

bin  q  "  0 

do  i  I ,  nuui  bins 

dischargc3(i)  ((heada  3(6)’''north  vel3(i))-(hcada  3(7)* 
I  east  vcl3(i)))*(bin  size/l00)*units*dt 


if(abs(discharg  '3(i))  .It.  10000)  then 
bin  q  discharge3(i)  +  bin  q 
endif 
enddo 

headb  3(1 1)  =  bin  q  +  bin  total 

headb  3(3)  "=  heada_3(7)  *  dt  ^  n  total 

headb  3(4)  =  heada  3(6)  *  dt  +  e  total 

headb  J(5)  -  (headb_3(3)**2  +  headb  3(4)**2)* *0.5 

headb_3(l)  =  ((heada  3(7)*dt)**2  +  (heada_3(6)*dt)**2) 

+**0.5  + 1  total 

headb_3(12)  -  topbttn3(l)  +  last  top 
headb  3(13)  =  topbtm3(2)  +  last  btm 

write(2,  200)  (time3(i),  i=l,  7),  ens_num3,  (heada_3(i),  i=l,  13), 
+(bdepth3(i),  i=l,  4),  (headb_3(i),  Hi,  19),  labcl4 
do  j-1,  num  bins 

write(2,  300)  depth3(j),  vel  jtiag3(j),  vel  dirp3(j),  east  vel3{j), 

I  north  vel3(j),  vert_vel3(j),  error  vel3(j),  db  13(j),  db_23(j), 

+  db_330,  db_43(j),  per  j'ood3(j),  dischargc3(j) 
enddo 
ie3  =  ic3  -  1 
ie4  ~  ic4  1 

print  *,  'ensemble  no.  ',ens  nurn3,'  has  been  fixed' 

1  total  =  headb  3(1) 
n  total  -  headb  3(3) 
ejotal  =  headb_3(4) 
bin  total  headb  3(1 1) 
last  top  -  headb  3(12) 
last  btm  -  headb  3 ( 1 3) 

85  continue 

dl  l_total  -  headb  1(1) 

dn  =  ri  total  -  headb_l(3) 

de  =  e  total  -  headb_l(4) 

dbin  q  =  bin  total  -  headb  1(11) 

dtop  q  =  headb_3(  1 2)  -  headb  1(12) 

dbtm  q  -■=  headb  3(13)  -  headb  1(13) 

headb  2(  1 )"- headb  2(  1 )  t  dl 

headb  2(3)  -  headb  2(3)  t  dn 

headb _2(4)  headb_2(4)  -t  de 

headb  j2(5)  =  (headb  2(3)**2  -i  headb  2(4)**2)**0.5 

headb  2(  11)-=  headb_2(  1 1 )  t  dbin  q 

headb_2(12)  ”  headb_2(12)  +  dtop  q 

headb  2(13)  =  headb  2(13)  +  dbtm  q 

timcl  =  timc2 
cns  num  1  =  ens  num2 
heada  l  -=  heada  2 
bdepthl  “  bdepth2 


hcadb  1  head!)  2 

depth  1  dcpth2 
vcl  inagl  vcl  niag2 

vH  dim  I  vcl  dini2 

cast  veil  east  vel2 

nortli  veil  iiorlli  vel2 
vert  veil  vert  vel2 

error  veil  error  vel2 
db  11  db  12 

db  21  db  22 

db  31  db  32 

db  41  db_42 

per  good  1  pei  jjood2 
discharge  1  dischargc2 
last  t  hcadb  1(2) 
ie2  0 

goto  35 

continue 


end 


subroutine  clcarscrccii 
do  k  1, 25 
writc(*,  '• ) 
onddo 
end 


A1M*KNI)IX  1) 


GI’SQ  Version  I.O:  User  Maniial 

INTRODUCnON 

riiis  manual  describes  vei  sion  1  of  a  program  designed  to  utilise  data  from  an  external 
DilVerential  (ilobal  Positioning  System  to  re-calculate  river  discharges  measured  using  an  RDi 
instinnienls  ADCP, 

The  source  code  for  the  eomputei  program  is  written  in  FORTRAN  and  is  listed  in 
Appendix  B  An  executable  version  of  the  code  is  contained  on  a  MS-DOS  fomialted  disk, 
labelled  "GPSQ  Ver  1 .0".  The  program  will  lun  on  most  IBM  compatible  P.C.s.  This  program 
may  be  Ireely  copied,  but  the  source  code  may  not  be  distributed  to  other  parties.  The  author 
aecc{)ts  no  responsibil^  V  or  liability  resulting  from  the  use  ol  the  program. 

I^STALLAT!ON 

To  install  CJPSQ  ver  1 .0  onto  your  hard  drive,  insert  the  diskette  provided  into  the  a:  drive 
of  your  machine  luid  type  a:\in.stall.  A  directory  named  GPSQ  will  then  be  created  on  your  hard 
drive,  the  program  files  will  then  be  copied  into  this  directory. 

RllNiNl^t:  mu  PROGRAM 

Run  the  progiam  by  typing  GPSQ  at  the  command  prompt  in  the  GPSQ  directory.  A  title 
screen  should  then  be  displayed,  press  HN  THR  to  continue.  The  program  then  runs  through  the 
following  stages: 

( I )  rhe  user  is  prompted  for  the  name  of  the  ADCP  and  Navigation  (GPS)  data  files, 
if  the  lilcs  at  e  not  in  the  current  directory  then  a  path  must  be  specified.  Sample 
data  files,  ADt’P.DA'f  and  GPS. DAT,  are  contained  on  the  installation  disk. 


(2)  The  user  is  asked  to  input  a  name  for  the  final  discharge  data  output  tile.  11’  no 
path  is  specified,  the  file  is  put  placed  in  the  current  directoi-y 

(3)  A  menu  is  displayed  on  the  screen,  giving  options  for  data  output  units,  choose 
from  either  English  units  (ft,  ft/s,  efs)  or  SI  units  (m,  m/s,  cumocs). 

(4)  A  second  menu  is  displayed,  asking  the  user  to  input  a  choice  for  the  method  to 
be  used  for  top  layer  velocity  estimation,  either  "power"  or  "constant". 

(5)  A  'depth  calculation  options'  menu  is  displayed,  either  ADCP  beam  depths  or 
FATHOMETER  data  can  be  chosen.  At  present  the  fathometer  depth  data  option 
is  not  available,  and  if  choosen  an  c'^^'^r  message  to  this  affect  will  appear. 

(6)  The  program  prompts  the  user  for  the  ADCP  beam  angle  (in  degrees).  If  this 
information  is  inknown,  it  can  be  determined  by  pressing  F9  in  the  RDI 
PLAYBACK  program. 

(7)  The  program  asks  for  the  local  magnetic  declination  in  degrees,  if  not  already 
corrected  for  in  the  ADCP  configuration  file. 

(8)  A  bottom  track  declination  is  then  displayed  on  the  screen.  The  bottom  track 
declination  is  equal  to  the  average  directional  dilfei  ence  between  the  bottom  track 
data  and  GPS  data.  Hit  ENI’ER  to  continue 

(9)  A  power-law  co-efficient  is  then  required  by  the  program,  a  value  of  0, 1667  is 
recommended  for  most  applications,  for  more  information  see  Chen,  (1991). 

( 1 0)  There  is  then  a  short  delay  whilst  the  program  processes  the  velocity  data.  A  menu 
is  then  displayed  asking  the  user  whether  or  not  the  require  the  near  shore 
discharge  calculation  option. 

(1  ly  The  final  results  of  the  discharge  calculation  are  then  displayed  on  the  screen  in 
teiTHs  of  total  discharge,  total  area,  average  velocity. 

( 1 2)  During  the  program  two  files  are  created,  a  velocity  data  file  and  discharge  data 


file.  The  format  of  these  files  are  described  in  tlic  main  body  of  this  report. 


APPENDIX  E 


ADC'PFIEL:  User  IVIanual 

INTRODUCTION 

I’his  manual  describes  a  program  designed  to  fill  in  missing  data  from  processed  ADCP 
data  files,  in  the  RDl  TRANSHCT  program  fomial 

1  he  source  code  for  the  computer  program  is  written  in  FORTRAN  and  is  listed  in 
Appendix  C,  An  executable  version  of  the  code  is  contained  on  a  MS-DOS  foimatted  disk, 
labelled  "ADC’PFILT".  The  program  will  run  on  most  IBM  compatible  P  C  s.  This  program  may 
be  fieely  copied,  but  the  source  code  may  not  be  distributed  to  other  parties.  The  author  accepts 
no  responsibility  or  liability  resulting  from  the  use  of  the  program. 

INSTALUAIION 

To  install  AIXIPFILL  onto  your  hard  drive,  insert  the  diskette  provided  into  the  a:  drive 
of  your  machine  and  type  a:\install.  A  directory  named  ADCPF'II  ,L  will  then  be  created  on  your 
hard  drive,  the  program  files  will  then  be  copied  into  this  directory. 

UlINNINC;  THE  PROCiRAM 

Run  the  program  by  typing  ADCPFII.L  at  the  command  prompt  in  the  ADCPFII.L 
dircctoiy.  A  title  screen  should  then  be  displayed,  press  EN  TTR  to  continue.  1  he  program  then 
runs  through  the  following  stages: 

(1)  The  user  is  prompted  for  the  name  of  the  processed  ADCP  data  files,  if  the  files 
arc  not  in  the  current  directory  then  a  path  must  be  specified. 

(2)  The  program  then  processes  the  ADCP  file,  checking  for  bad  ensembles,  i'hc 
numbers  of  any  bad  ensembles  are  displayed  on  the  screen  at  this  time. 


A  new  ADCP  data  lile,  CLEAN. OUT,  with  all  bad  or  missing  data  segments  filled 
in  is  produced  in  the  current  working  directory. 


